The Heavy Quark Spin Symmetry Partners of the X(3872) 
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We explore the consequences of heavy quark spin symmetry for the charmed meson-antimeson 
system in a contact-range (or pionless) effective field theory. As a trivial consequence, we theorize 
the existence of a heavy quark spin symmetry partner of the X(3872), with J^^ = 2+"*", which we 
call X(4012) in reference to its predicted mass. If we additionally assume that the X(3915) is a 0^"'' 
heavy spin symmetry partner of the X(3872), we end up predicting a total of six _d(*)_d'*^ molecular 
states. We also discuss the error induced by higher order effects such as finite heavy quark mass 
corrections, pion exchanges and coupled channels, allowing us to estimate the expected theoretical 
uncertainties in the position of these new states. 

PACS numbers: 03.65.Ge, 13. 75. Lb, 14.40.Lb, 14.40.Pq, 14.40.Rt 



I. INTRODUCTION 

The discovery of the X(3872) resonance [l] might 
have confirmed a well-known theoretical expectation of 
hadronic physics, heavy meson molecules The 
X(3872), with a mass of mx = 3871. 57 ± 0.25 MeV 0, is 
extremely close to the £)°*D° threshold {rrijjo + mjj.o = 
3871.73 ±0.21 MeV), a feature suggesting that its nature 
is mostly molecular, where the role of other more com- 
pact components (e.g. tetraquark, cc) will be less im- 
portant in comparison. The J^*^ quantum numbers of 
the X(3872) are either 1++ or 2"+ ^8.-10,]. with a slight 
experimental preference for the second set of quantum 
numbers However, only the 1+"'" assignment is com- 
patible with the molecular interpretation. If the 2 ^ 
possibility turns out to be the correct one, the X(3872) 
would be much more exotic than expected. In this re- 
gard, the analysis of Hanhart et al. [ll[ suggests that 
the experimental data, though too scarce to draw defini- 
tive conclusions, might indeed be more compatible with 
1++ after all (see, however, Ref. [l3| for a different opin- 
ion). Meanwhile, in the bottom sector, the recent dis- 
covery of the Zf,(10610) and Zb(10650) isovector states 
by Belle [H, [ll| also provides two new strong candidates 
for molecular states, as the Zi,'s lie close to the B*B and 
B*B* thresholds respectively. 

Actually, not all the heavy meson molecules are neces- 
sarily as shallow as the X(3872) or the two Zi, resonances. 
Several theoretical works have suggested the molecular 
interpretation of other XYZ states: the X(3915) [lH| 
and the_r(4140) [ll| have been identified as a D*D* 
and DID* bound state respectively by different theoret- 
ical approaches [i3-[ll- The y(4260) ^ has eveii been 
proposed to have the three body structure J/'^KK [2l[. 
Moreover other states have been predicted, especially in 
the isoscalar bottom sector [22 - l24| , but have not yet been 
experimentally confirmed or discarded. 
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Heavy meson molecules are a natural thing to ex- 
pect on the basis of the similarity between the meson- 
meson interaction and the nuclear force that binds the 
deuteron 0, [ij. From this analogy, we expect that 
the effective field theory (EFT) formulation of nuclear 
physics [25-29] will also represent a constructive ap- 
proach to the description of heavy meson systems at 
low-energies. As in the nucleon-nucleon system, the 
low energy interaction between a pair of heavy mesons 
is mediated by pion exchanges, which in turn are con- 
strained by chiral symmetry. In contrast, the nature of 
the short range interaction remains unknown, but we can 
parametrize it in terms of contact-range operators be- 
tween the nucleon / heavy meson fields. However, in 
the case of heavy meson molecules there is a particularly 
simplifying feature: pion exchanges are weaker than in 
the nuclear case, owing to the smaller light quark con- 
tent of the heavy mesons in comparison to the nucle- 
ons. This means that pions are amenable to a perturba- 
tive treatment in a larger range of energies than in the 
nucleon-nucleon system. From this it is expected that the 
EFT description of heavy meson molecules will simplify 
at lowest order to a contact range theory, at least for a 
certain binding energy window 1301. A nice illustration of 
this idea is provided by X-EFT [3l|, which considers the 
low energy description of the X(3872) state as a D^D^* 
molecule. 

On a different level, the presence of a charm 
quark/antiquark in the heavy meson and antimeson con- 
forming the X(3872) dictates that heavy quark spin sym- 
metry (HQSS) [33 - l35| is relevant for this system. In the 
context of the EFT description of heavy meson molecules, 
HQSS constrains the form of the contact range oper- 
ators of the theory in a very specific way ^^\. What 
this means is that there should be HQSS partners of the 
X(3872), in analogy with the theorized HQSS partners of 
the .^{,(10610) and Zh(10650) that have been already dis- 
cussed in Rcfs. [37, 38] . The purpose of this paper is to in- 
vestigate the HQSS structure of charm meson-antimeson 
(D, D* and D, D*) molecules, with the intention of iden- 



2 



tifying the possible HQSS partners of the X(3872). For 
that we will assunie that certain XYZ states are molec- 
ular, in particular the X(3915) [l^. This identification 
will span a total of six states, some of which have also 
been predicted in other schemes. 

The article is structured as follows: in Sect. [II] we 
present the EFT description of heavy meson molecules 
that we advocate, which is in turn inspired on the ideas 
of Ref. lOl. We explore the consequences of this EFT 
in Sect. IIIIl where we deduce from HQSS the existence 
of the 2++ D*D* partner of the X(3872), and of other 
four additional states if we identify the X(3915) as a 0++ 
D*D* molecule. In Sect. lIVl we probe the robustness of 
the previous results by considering the effect of sublead- 
ing order contributions, namely finite heavy quark mass 
corrections, pion exchanges and particle coupled channel 
effects. In Sect. El we briefly discuss some HQSS con- 
straints on the total decay widths of the states found in 
this work. In Sect. IVll we discuss the significance of the 
results. Finally, the technical details of the manuscript 
can be consulted in Appendices El IE) and O 



II. THE EFT DESCRIPTION AT LOWEST 
ORDER 

In this section we review the EFT framework that 
we use for the description of heavy meson-antimeson 
molecules at lowest order. The presentation is simple and 
schematic, centered in the conceptual issues rather than 
the specific details, which can be consulted in Ref. (30| . 
The EFT we are advocating is in fact identical to the 
one presented in Ref. ^] for the isovector bottom sec- 
tor. Following the findings of Ref. [13], we assume pion 
exchanges and particle coupled channel effects to be a 
subleading correction entering at next-to-leading order 
(NLO) and next-to-next-to-leading order (N^LO) respec- 
tively. Nevertheless, we will perform explicit calculations 
to test these assumptions. 



A. Overview of the EFT Formalism 

The EFT approach provides the possibility of con- 
structing generic and systematic descriptions of arbitrary 
low energy processes. They are particularly useful when 
the system we are interested in cannot be easily explained 
in terms of a more fundamental description at higher en- 
ergies. The EFT idea is simple: we identify the fields 
and symmetries that are relevant at low energies and 
construct all possible interactions compatible with them. 
Even though the number of interactions is infinite, they 
can be classified according to their importance at low en- 
ergies by means of power counting, the ordering principle 
of EFT. If Q is the soft (low energy) scale of the system 
we are describing and Aq the hard (high energy) scale, 
power counting allows us to express any physical quan- 
tity as a power series in terms of the small parameter 



Xq = Q/Aq. 

For illustrating this idea let us consider a physical 
quantity A that we want to compute in the EFT frame- 
work. This quantity receives in principle contributions 
from all the relevant diagrams involving the low energy 
fields and compatible with the low energy symmetries: 



A(g,Ao) = ^A(^)(g,Ao), 



(1) 



However, the different diagrams have different scaling 
properties that can be used for ordering the sum above. 
For example, we have the canonical dimension of A, 
which is defined as 



A(^) (A Q, A Ao) = A''^ A(^) (Q, Ao) 



(2) 



and is the same for all the EFT contributions to A. But 
the interesting scaling property is power counting, which 
refers to the behaviour under a transformation of the type 
XQ 



A(^)(Ag,Ao) = A''- A(^)(Q,Ao). 



(3) 



where i^o is the order of the contribution D, which is 
bounded from below (i.e. > i^a)- The sum of diagrams 
above can be reorganized as an expansion in terms of 
increasing scaling dimension: 



A(Q,Ao)= J2 A^'^HQ.^o), 



(4) 
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where, for simplicity, we have get rid of the D subscripts 
and superscripts. For each order there is only a finite 
number of diagrams that contributes to the quantity A. 
Combining the scaling laws of Eqs. ([2]) and ([3]) we obtain 
a well-defined power series for A 

A(g,Ao)=A^- y: 



(5) 



with A*-"^ a dimensionless quantity that we expect to be 
of the order of unity (i.e. Q^). Notice that A^^^ does not 
depend on the hard scale Aq and is related to A'^^^ via 
Eqs. ([2]) and ([3]). In the formula above, Q' is an auxiliary 
soft scale we use to express as a function of a dimen- 
sionless ratio. Provided there is a clear scale separation 
in the system, that is, Aq > Q, the power series above 
will be convergent. Not only that, if we consider only 
contributions from diagrams with v < fmax, the error of 
the EFT calculation will be Xq"""'"''^. In this work we will 
only perform calculations at the lowest order and we ex- 
pect a relative error of the order of xq in the calculations 
to follow. 

If we are interested in the low energy description of 
heavy meson-antimeson bound states, the relevant phys- 
ical object we want to expand is the (non-relativistic) 
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potential between the heavy meson and antimeson: 



(6) 



The expansion starts at order vq > —1, where xq is the 
ratio of the soft and hard scales of the system. The low 
energy degrees of freedom we consider are the heavy me- 
son and antimeson fields and the pion field. The pion- 
meson vertices are constrained by chiral symmetry and 
the corresponding Feynman rules can be derived from 
heavy hadron chiral pertm-bation theory (HHChPT) [s^ . 
In turn HQSS generates strong constraints on the form 
of the heavy meson-antimeson interactions [36] . This 
means that the EFT potential includes two kind of con- 
tributions: contact range interactions, i.e. four meson 
vertices, and pion exchanges. The set of soft scales Q 
includes in principle the pion mass m^r and the center-of- 
mass momenta p and p' of the meson and antimeson. The 
hard scale Aq can represent the momentum scale at which 
we expect the low energy symmetries to break down, i.e. 
the chiral symmetry breaking scale A^^ — A-Kf.,^ ~ 1 GeV 
(with the pion decay constant) for chiral symmetry 
and the heavy quark mass mg for HQSS, but it can also 
stand for the momentum scale at which the composite 
structure of the heavy mesons starts to be resolved. 

At lowest (or leading) order {Q^) the heavy meson- 
antimeson EFT potential is local and only depends on 
the momentum exchanged by the meson and anti-meson 



(P|1/(")|P')=^^°'(P-P'), 



where we have 



q ■ 



(7) 



(8) 



As can be seen, the potential is the sum of a contact and 
a finite range contribution. The contact range operator 
Co is a free parameter of the theory. The finite range 
contribution is the well-known one pion exchange (OPE) 
potential, where g ~ 0.6 is the axial coupling between 
the heavy meson and the pion (we have particularized its 
value for the charmed meson case, see Refs. [i^ |4l| for a 
determination), ~ 132 MeV the pion decay constant, 
and q the momentum exchanged by the heavy meson and 
antimeson. The sign 77 and the polarization operators a 
and b depend on whether the initial/final states is PP, 
PP*, P*P or P*P*. A more detailed account can be found 
in Appendices [XI IB] and O where the leading order (LO) 
potential is derived. 

A problem with the EFT potential above is its be- 
haviour for larges values of the exchanged momentum, 
|g I 3> TOjr. In this limit the LO potential tends to a 
constant value. At higher orders in the EFT expansion 
the problem worsens and the potential actually diverges. 
This feature can be easily deduced from power counting 



which admits a solution of the type V^'^\q) oc \q\'^ for 
|(f I ^ TOtt. Of course, if we are only considering tree level 
amplitudes involving no heavy meson-antimeson loops, 
then there is no conceptual problem with the previous 
divergences: the EFT potential is only expected to make 
sense at low energies. However, if we iterate the potential 
in the Schodinger or Lippmann-Schwinger equation (a 
necessary step for the description of bound states) the 
divergences will require renormalization. 

In non-relativistic EFTs the renormalization process 
is straightforward. First, we begin by regularizing the 
potential: 

{p\VA\p')=f{^){p\V\p')fi^), (10) 

where Va is the regularized potential that we will employ 
in actual calculations, V is the unregularized (i.e. the 
original) potential, f{x) a regulator function and A an 
ultraviolet cut-off. In the calculations to follow, we will 
use a gaussian regulator of the type f{x) = . At this 
point physical predictions can still depend strongly on 
the value of the cut-off, a problem we must solve. Thus 
there is a second step in the renormalization process: we 
allow the contact range operators to depend on the cut- 
off. At lowest order this implies that (7^°) = C(°)(A). If 
we have included all the counterterms required by power 
counting, they will be able to absorb all the divergences 
of the theory. The calculations will still contain a residual 
cut-off dependence, but this does not represent a prob- 
lem: its size is expected to be a higher order effect, at 
least for a judicious choice of the cut-off. Even though 
there is not a well-established criterion for choosing the 
cut-off, calculations in nuclear EFT suggest that the opti- 
mal value of the cut-off should never be much larger than 
the hard scale Ag [12] . Here, for avoiding the breakdown 
of the low energy symmetries in loops, we advocate the 
slightly more stringent condition A < Aq . We employ the 
actual values A — 0.5 GeV and 1 GeV, for which we have 
checked that the results do not change considerably. This 
indicates that renormalization has been correctly imple- 
mented. 



B. Power Counting and Bound States 

The description of bound states requires the iteration 
of the EFT potential in the bound state equation 



(11) 



where \^ b) is the wave function, Go{E) — l/{E~Ho) the 
resolvent operator and V the potential. If we require the 
bound state equation to be compatible with the power 
counting of the EFT potential, then successive iterations 
of the GqV combination must be of the same order: 



T/M(Ag) = A'^t/('')(Q), 



(9) 



0{V) =0{VGoV). 



(12) 
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By taking into account that the Go operator scales as Q 
in loops we can see that only the order contribu- 
tion to the potential should be iterated. Thus the pres- 
ence of shallow bound states or large scattering lengths in 
a two-body system requires the non-perturbative treat- 
ment of a piece of the effective potential. There ex- 
ists a problem then: the EFT potential we obtain from 
HHChPT starts at order and is therefore incompatible 
with the EFT description of a low energy bound state. 

The solution is to redefine power counting by pro- 
moting the Cq contact range operator from order to 
|43l - l47| . This is equivalent to assume that the Co 
operator is contaminated by a low energy scale. We can 
confirm this assumption a posteriori by solving the bound 
state equation with the Co operator alone: if we regular- 
ize the EFT potential with a cut-off A, and set the value 
of Co (A) to reproduce the position of the bound state, 
we obtain the generic result 



1 



Co (A) 



27r 



7B A 

TT 



(13) 



where /i is the reduced mass of the two-body system and 
7b = \/^2fj,EB the wave number (with Eb the bound 
state energy). Of course, the exact form of the relation 
depends on the specific regularization scheme employed 
- the expression above corresponds to power divergence 
subtraction [i^, |3] ~ but will be in the line of the pre- 
vious form. From the equation above it is evident that 
we can only have Co ^ as long as the cut-off A 

scales like 0{Q). Note that this is consistent with our 
requirement that all involved momenta should be soft 
and smaller than the charm quark mass to make sense of 
HQSS. 



C. The EFT Potential at Lowest Order 

As we have seen, the iS-wave LO interaction {0{Q^^)) 
between a heavy meson (H = P,P*) and antimeson 
(H = P,P*) only contains contact operators (i.e. four 
heavy meson vertices). The contact range interac- 
tions arc in turn constrained by HQSS and depend 
on the particle channel under consideration (HH — 
DD,DD*/D*D,DD), the total isospin / and the value 
of the J^'^ quantum numbers. In particular HQSS lim- 
its the number of independent LO counter-terms to two 
per isospin channel at LO [36j . We will only consider the 
isoscalar / = channels, to which the AT (3872) and the 
^(3915) belong, in the isospin symmetric limit. 



^ This can be trivially checked by considering the rescaling trans- 
formation 

where the energy rescales as as we are considering a non- 
relativistic system. 



The HQSS contact range interaction can mix different 
particle channels with the same J^*^ quantum numbers. 
Therefore, for writing the LO potential we considering 
the set of particle coupled channel basis 

6(0++) = {|PP),|P*P*(0))} , (14) 
6(1+-) = {\PP*) + \P*P)) , |P*F*(1))| ,(15) 

S(l++) = |-i=(|PP*)-|P*P))| , (16) 
B{2++) = {\P*P*{2))} , (17) 



where the number in parenthesis in the |P*P*(S')) states 
is the total intrinsic spin S to which the vector meson- 
antimeson system couples. In this basis, the EFT poten- 
tial is independent of momentum and reads (see Appen- 
dices EE] and [C|) 

^^°teo++) = f ) , (18) 



, V (19) 

Oa — Cob / 



F^^(g,l++) = Coa+Cob 



Ob : 



(20) 
(21) 



where Coa and Cob are the two independent counter- 
terms that we expect from HQSS for a given isospin 
sector. This contact potential V behaves as 0{Q^^), 
but for avoiding confusions with the notation we have 
renamed to V^'-', as the negative exponent could 

be confused by the inverse of the potential. As can be 
appreciated, the 1++ and 2++ cases are uncoupled and 
their potential is identical, a strong hint that we should 
expect a 2++ D*D* partner of the X(3872). 

At this point we notice that the heavy pseudoscalar 
and vector mesons P and P* are only degenerate in the 
heavy quark limit mg — > oo. For finite rng there is a 
mass splitting between the heavy mesons 



Mf 



Mf 



Ar 



(22) 



that scales as l/mg. As a consequence of this gap, the 
two HH thresholds in the 0++ and 1^ coupled channel 
happen at different energies. If we are interested in low- 
lying bound states, the energy difference between the two 
thresholds may actually be considerably larger than the 
binding energy of the state. Within the EFT framework 
this means that we may very well be entitled to ignore the 
coupled channel effects. The momentum scale associated 
with the coupled channels is 



Ac(0++) = .\/2/i(2AQ), 

Ac(l+-) = V2jrA^, 



(23) 
(24) 



for the 0++ and 1^ cases respectively, where /i is the 
reduced mass of the HH heavy meson system and Ag 
is the energy gap. In the particular case of the charm 
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mesons, direct evaluation yields Ac(o+^ 



A 



c(i+-) 



520 MeV. If we have a 



++ 



- 750 MeV and 
DD (1+^ D*D) 



bound state at threshold, the corresponding 0++ D*D* 
(1^' D*D*) component will have a wave number at least 
of the order of the hard scale Aq ~ 0.5 — 1 GeV. Thus 
there is no problem in ignoring the coupled channel struc- 
ture and treating the two particle channels in the 0"*"^ 
and 1"' cases as independent. From a more formal EFT 
viewpoint what we are doing is to consider the coupled 
channel momentum scale as Ac ^ 0{Q^), from which we 
expect the Gq operator involved in the particle mixing to 
scale like (see for instance Ref. [131). This translates 
into a suppression of particle coupled channel effects by 
two orders in the chiral expansion: if we count the Cob 
operator as Q^^, then particle coupled channels do not 
enter until order Q, that is, at least one order beyond 
pion exchanges. 

Taking into account that we can distinguish between 
different particle channels in the EFT, the LO potentials 
for the 0++ and 1^ cases finally simplify to 





) = Coa 


5 


(25) 




) = Coa 


— 2 Cob , 


(26) 


T/LO _ (ff 

^p*p/pp* 


) = Coa 


— Cob , 


(27) 




) = Coa 


— Cob , 


(28) 



where we have added the particle channel as a sub- 
script for distinguishing states with the same J^*^ quan- 
tum numbers. Even though the particle coupled channel 
structure of the 0++ and 1"' molecules can be ignored 
in the LO description of these states, we expect the O"*"^ 
and 1"' p*p* states to have a strong tendency to decay 
to PP and PP* respectively. On the contrary, the 2++ 
P*P* will have a smaller partial decay width to a heavy 
meson-antimeson pair, as this process does not happen 
via a counterterm, and therefore will be suppressed 
by one order in the EFT expansion. 



D. Bound States at Lowest Order 

Finally, for completeness, we briefly discuss the solu- 
tion of the bound state equation for the contact range 
potentials that appear in the LO description of heavy 
molecular states. Even though well-known, it will be of 
help for the calculations in the next section. We solve 
the bound state equation with the regularized potential 
Va that reads 



(p|^A°|p")=/(?)Co(A)/(?), 



^A' 



-A^ 



(29) 



where f{x) is the regulator, and Co the appropriate coun- 
terterm for each J^*^ / particle channel combination. For 
this potential, the wave function admits the ansatz 



f(P^ 
p + 7 A 



(30) 



where A/" is a normalization constant, fi the reduced mass 
of the two body system and 7 = ^—2^Eb the wave num- 
ber of the state, with Eb the binding energy. Direct sub- 
stitution into the bound state equation (Eq. (dJ)) yields 



Co (A) 



(fig 



g 2^ 



(31) 



This is the eigenvalue equation of the contact range the- 
ory. For the HH system there are six of these equations 
corresponding to the six possible S-wave states: in each 
case, we simply need to particularize the values of Co and 
fi as appropriate. For further details we refer the reader 
to Ref. [2J|, where the equations above were derived 
and discussed in the context of heavy meson-antimeson 
molecules. 



III. THE PARTNERS OF THE X(3872) 

We start by considering the 1++ DD* state, the 
AT (3872), within the EFT formalism described in the pre- 
vious section. From HQSS we expect that the heavy 
meson-antimeson interaction in the 1^^ DD* and 2+''' 
D*D* channels will be identical. In terms of the EFT 
potential we have 



V 



LOQ++ 



(32) 



as deduced from Eqs. (|20ll2ip . From this, we automat- 
ically anticipate the existence of a isoscalar 2++ D*D* 
bound state with a binding energy similar to that of the 
A:(3872). 

To pinpoint the exact location of the 2++ partner of 
the A" (3872) we begin by determining the counterterm 
combination Coq + C06. We will use a gaussian regu- 

2 

lator f{x) = e^^ and the cut-off values A = 0.5 GeV 
and 1 GeV. As we are working in the isospin symmetric 
limit 0, we can consider the X(3872) to be a D*D bound 
state with a binding energy of Bx{l'^~^) — 4.2 MeV, from 
which we obtain Coa + Cob = -1.94 fm^ (-0.79 fm^) for 
A = 0.5 GeV (1 GeV). Now we can predict that the mass 
of the 2++ state lies in the vicinity of 4012 MeV (a value 
rather independent of the cut-off), corresponding to a 
binding energy of 5 MeV. We call this state the Ar(4012), 
and stress that this prediction is independent of any as- 
sumptions about the molecular nature of any other XYZ 
states, relying on HQSS alone. We mention in passing 
that isospin breaking effects, even though crucial for un- 
derstanding certain decay properties of the Ar(3872) (see, 
for example, Refs. 0,1431), ^.ot an appreciable effect 
in the spectroscopy problem 0. 



In what follows, 



we use the isospin averaged masses tuq 
2008.6 MeV. 



1867.2 MeV and m^j 
^ If we take into account isospin breaking within the formalism of 
Ref. [49I ]. we find that the position of the Jsr(4012) moves by about 
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HH 




Vc 


E {A^ 0.5 GeV) 


E {A^l GeV) 


Exp [7] 


0++ 


DD 


'So 


Coa 


3706 ± 10 


3712l:J-^ 


- 


1++ 


D*D 


'Si 


Coa + Cob 


Input 


Input 


3872 


1+- 


D*D 


^Si 


Coa — Cob 


3814 ± 17 


3819^27 




0++ 


D*D* 


'So 


Coa — 2 Cob 


Input 


Input 


3917 


1+- 


D*D* 


'Si 


Coa — Cob 


3953 ± 17 


3956^2! 


3942 


2++ 


D*D* 


'S2 


Coa + Cob 


4012 ± 3 


401219 





TABLE I. Predicted masses (in MeV) of the X(3872) HQSS partners for two different values of the gaussian cutoff. We use 
as input 3871.6 MeV and 3917.4 MeV for the X(3872) and X(3915) masses, respectively. From this we find that the value of 
the LO couplings are Coa = -3.53 fm^ and Cot = 1-59 fm^ for A = 0.5 GeV and Coa = -1.06 fm^ and Cot = 0.27 fm^ for 
A = 1 GeV. Errors in our predicted masses are obtained by varying the strength of the contact interaction in each channel by 
=F15%, which corresponds to the expected violations of HQSS for the charm quark mass. 



For predicting states beyond the X(4012), we have to 
identify a particular XYZ state as a further molecular 
partner of the X(3872). In this way we will be able 
to determine the two contact range couplings (Coa and 
Cob) and obtain the full spectrum of molecular states. 
Two interesting candidates are the X(3915) [l^ and 
the X(3940) [sol . [Slj . from which the first is the most 
promising. The X(3915) has been theorized to be a 0++ 
or 2++ D*D* molecule in Refs. [13, Even though 

the work of Refs. [13, [l^l cannot discriminate between 
the 0++ or 2++ quantum numbers 0, HQSS suggests 
that the most probable J^*^ value is 0++ instead of 
2++, as the later would imply a remarkable violation 
of HQSS. On the contrary, the accommodation of the 
X(3940) 13 Hi within the HQSS pattern of molecu- 
lar states faces a problem. The X(3940) decays strongly 
to DD*, a feature compatible with the expected prop- 
erties for a l^ D*D* axial state. However, the pro- 
duction mechanism for the X(3940) is more compati- 
ble with a positive C-parity state than with a nega- 
tive C-parity one; this state is produced in the reac- 
tion e+e^ J/\I'X(3940), most probably via an inter- 
mediate virtual photon (e+e~ — > 7* — >• J/\l/X), sug- 
gesting that the C-parity is positive. Nonetheless this 
is not a definitive conclusion, and it may happen that 
the X(3940) resonance is bein g pr oduced via two virtual 
photons, see for example Refs. j52l.[53| for a case in which 
this process is less suppressed than naively expected. 

At this point we notice that the identification of the 
X(3915) as a D*D* molecular state, though more promis- 
ing than that of the X(3940), is not free of problems 
either. In particular, the binding energy of the corre- 
sponding heavy vector meson-antimeson system is of the 
order of ~ 100 MeV, which translates into a wave num- 



~ 1 MeV, to 4013 MeV. In this calculation we have ignored the 
1 = 1 counterterm in the J^^ = 1++ channel, which (in absence 
of additional information) is considered to be of order and 
hence subleading. 

Branz et al. llSlI notice that the (little) known decay properties 
of the X(3915) are compatible with both assignments 



ber of 450 MeV. This means that the X(3915) lies not 
too far away from the limits of what can be described 
within the EFT. Its wave number indicates that a de- 
scription in terms of mesons alone may be incomplete 
and that the explicit inclusion of shorter range compo- 
nents (e.g. tetraquark or charmonium-like) may be neces- 
sary [3. However, the EFT framework is very helpful and 
convenient in this regard. Abusing the limits of the EFT 
translates into a noticeable cut-off dependence and a lack 
of convergence of the EFT expansion, that is, subleading 
order corrections will be able to completely alter the LO 
results. As we will explain in the next paragraph, the cut- 
off dependence is numerically small, and as we will check 
in the next section, the subleading order corrections are 
moderate, but nonetheless under control. All this indi- 
cates that the X(3915) is probably more amenable to an 
EFT treatment than naively expected. 

From the assumption that the X(3915) is a 0++ 
molecule we can determine the counterterm combination 
Coa — 2 Cob and consequently the location of the six pos- 
sible HQSS partners of the X(3872). The masses of the 
molecular states resulting from the previous identifica- 
tion can be consulted in Table HI We obtain a 0++ DD, 
l+~ DD* and l+~ D*D* state which we call X(3710), 
X(3820) and X(3955). The errors in Table |T] refer to 
uncertainties owing to violations of HQSS in the charm 
sector (see Section HVAI for a detailed explanation). The 
location of these three hadronic molecules is rather inde- 
pendent of the cut-off. If (instead of a gaussian regula- 
tor) we use a sharp cut-off, there are small variations of 
about ~ 1 MeV in the location of the states. Curiously, 
the X(3955) state we obtain is not far away from the 
aforementioned ^(3940) molecular candidate. A pos- 
sible identification is suggestive but contingent on the 
eventual determination of the quantum numbers (espe- 
cially the C-parity) of the X{3M0). 



This is based on the observation that the mean quadratic sepa- 
ration of the mesons is ^ (r^) = 0.5 — 0.8 fm depending on the 
cut-off. 
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E {A = 0.5 GeV) 


E {A^l GeV) 


Exp [7] 


0++ 


DD 


'So 


3708 


3720 


- 


1++ 


0*0 




Input 


Input 


3872 


1+- 


D*D 




3816 


3823 




0++ 


D*D' 




Input 


Input 


3917 


1+- 


D*D* 




3954 


3958 


3942 


2++ 


0*0* 


D2-^ S2- D2-^G2 


4015 


4014 





TABLE II. Predicted masses (in MeV) of the X(3872) HQSS partners when the OPE potential is included. We display results 
for two different values of the eraussian cutoff. Now, we find Coa = -3.46 fm^ and Cob = 1.98 fm^ and Coa = -0.98 fm^ and 
Cob = 0.69 fm^, for A = 0.5 and 1 GeV, respectively. 



IV. SUBLEADING ORDER CORRECTIONS 

In this section we explore the impact of the main sub- 
leading order contributions. For this we must take into 
account the existence of two different, unrelated expan- 
sions in the EFT formalism we are proposing. The first 
is the expansion in terms of inverse powers of the heavy 
quark mass and the second the standard power count- 
ing expansion. In the previous section we have assumed 
exact, instead of approximate, HQSS. The existence of 
1 /mg deviations from the heavy quark limit implies that 
the location of the molecular partners of the AT (3872) 
will be subjected to uncertainties. Apart from this, the 
power counting expansion indicates that we should take 
into account two important subleading order corrections 
appearing at order Q'^ and respectively: the OPE po- 
tential and the particle coupled channel effects. As we 
will see, the induced shifts in the position of the molecu- 
lar states from these subleading interactions will in gen- 
eral agree with (and in the case of the OPE potential be 
smaller than) the a priori EFT expectations, thus con- 
firming the robustness of the molecular spectrum we have 
deduced so far. 



A. The l/mq Corrections 

The LO potentials for the and 2++ heavy meson 
molecules - the X(3872) and the theorized X(4012) ~ 
are only identical in the heavy quark limit. Thus the 
existence of the X(4012) may be contingent on the size 
of HQSS violations stemming from the finite charm quark 
mass. In general, we expect the heavy meson-antimeson 
potentials of Eqs. (I20j|2ip and (|25ll28p to deviate from the 
heavy quark limit by a quantity of the order of 

= + (33) 

where rric is the charm quark mass (~ 1.5 GeV) and 
Aqcd ^ 200 MeV, translating into an expected 15% vio- 
lation of HQSS for the LO contact range potentials. The 
exceptions are the potentials of the 1++ DD* and 0++ 



D*D* channels, as they are fixed to reproduce the posi- 
tion of the X(3872) and X(3915) states. 

Assuming this 15% uncertainty in the contact range in- 
teractions, we obtain the error bars of Table ID As can be 
seen, the prediction of the 2++ partner of the ^(3872) 
is robust with respect to this theoretical error source. 
Other molecular states show moderate uncertainties of 
the order of 10 — 20 MeV in the binding energies. As 
we will see in the following subsections, these uncertain- 
ties are a bit smaller (but yet of the same order) as the 
subleading order corrections coming from pion exchanges 
and the particle coupled channels. 

B. One Pion Exchange 

A possible issue with the present EFT treatment of 
the X(3915) as a D*D* molecule is whether the OPE 
potential is really perturbative in this case. With a 
wave number of 7 ~ 450 MeV, the wave function of 
the Ar(3915) probes the intermediate distances where the 
tensor component of the OPE potential is stronger than 
the central one. In this regard and according to the ar- 
guments in Ref. we expect tensor OPE to become 
non-perturbative at a critical center-of-mass momentum 
of fccrit ^ 420 MeV for the 0++ D*D* channel 0. Con- 
sequently the Ar(3915) lies at the edge of the domain of 
validity of the EFT description we are using, in which 
pion exchanges are perturbative (i.e. small), and may 
require a more sophisticated EFT with non-perturbative 
pions. If we were predicting the binding energy of this 
molecular state, the previous observation would translate 
into a large uncertainty in the calculations of the order of 
100%. The cause of this uncertainty will be high energy 
fluctuations in the meson-antimeson loops generated by 
the pion exchanges. However, as the binding energy is 
used as the input of the calculation, the uncertainty can 
manifest either in the value of the Coa and Cot countert- 
erms, thus subjecting the predictions of Table |T] to large 



As we will comment later, this number is subjected to large un- 
certainties. 
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uncertainties, or in a failure of the theory to converge 
once subleading order corrections are included. 

Yet there are two mitigating circumstances that may 
increase the expected breakdown scale of the present 
EFT calculation, thus turning the predictions much more 
reliable. On the one side, the critical momenta above 
which the tensor OPE is no longer perturbative [3(?| are 
subjected to considerable uncertainties stemming from 
the far-from-perfect separation of scales in heavy meson 
molecules. In particular the value of the critical momenta 
may be up to 50% larger than expected, in which case 
the X(3915) will lie well within the range of applicability 
of the EFT with perturbative pions. On the other, we 
are limiting ourselves to the spectroscopy problem. It is 
worth noticing in this respect that the spin-spin structure 
of the Cob contact operator is very similar to the one we 
obtain from pseudoscalar meson exchange for the non- 
tensor piece of the interaction. We thus expect the Co& 
operator to be able to partially absorb the shift in the 
binding energy generated by the pion exchanges. How- 
ever, the effect is restricted to the central component of 
OPE, which is much better behaved than the tensor one 
at short distances. Of course, the definitive test is to re- 
calculate the position of the predicted molecular states 
after the non-perturbative inclusion of the OPE poten- 
tial. 

The distinctive feature of the OPE potential is that it 
can mix different partial waves. For the noii-pcrturbative 
calculation we follow the formalism of Ref . [24| , where we 
considered the partial wave decomposition of tensor OPE 
in detail for the PP* particle channel. The extension to 
the P*P* case is trivial and only entails a change in the 
partial waves that are actually coupled, which can be con- 
sulted in Table HT] For determining the value of the Coa 
and Cob counterterms we fix the binding energies of the 
X(3872) and X(3915) respectively. This procedure gen- 
erates again a total of six molecular states, four of them 
predictions, as can be seen in Table |IT1 We can appreci- 
ate that the binding energies of the states are relatively 
stable with respect to the iteration of the OPE poten- 
tial. The most affected state is the 0++ DD molecule, 
the X(3710) state, for which the binding energy can in- 
crease by almost 10 MeV for A = 1 GeV with respect 
to the pionless case. Curiously, for this state there is 
no OPE contribution to the finite range potential: two 
pseudoscalar objects cannot exchange a single pion, so 
the first non-trivial contribution to the finite range po- 
tential comes from two pion exchange, which we have 
not considered here. The dynamics of this state is solely 
controlled by the Coa counterterm. In this regard, when 
the contact range potential is adjusted to reproduce the 
binding energies of the X(3872) and X{3915) states in 
presence of the OPE potential, the change in the value 
of the Coa coupling is not counterbalanced by a pion ex- 
change contribution in the O"*"^ DD case. Hence, we find 
a larger shift in the location of this molecular state in 
comparison to the others. 

Indeed we find relatively small shifts in the energy of 



states other than the X(3710) (see Tables U and H!]) . The 
aforementioned observation seem to confirm the suspi- 
cions about the role of the Cob operator which could ef- 
fectively accommodate most of the OPE effects. Equiv- 
alently, the inclusion of the OPE potential produces a 
larger change in the Cob coupling than in the Coa one, 
and actually we find that the bulk of the change in Cob 
is given by the size of the contact piece of OPE 

9^ 

Cob ^ Cob — -TTn ' (34) 

where g'^/2f^ ~ 0.4fm^, see Appendix ICl for details 

The natural conclusion of the previous calculations 
is that OPE is indeed perturbative in heavy meson 
molecules. There is however a caveat: if the spectroscopy 
problem is rather insensitive to the OPE potential, how 
can we appreciate whether OPE is perturbative or not in 
the molecular states? A partial answer is that Cob is only 
expected to absorb effects coming from central OPE, but 
not from tensor OPE, which is the problematic piece. 
But there is another answer, that lies in the observa- 
tion that non-perturbative OPE would entail a signifi- 
cant change in the power counting of the contact range 
operators. As happens in nuclear EFT, the LO countert- 
erms that are able to renormalize the scattering ampli- 
tude with perturbative pions p3l |44| are not enough to 
renormalize the corresponding non-perturbative formu- 
lation (53 |. In particular, for heavy meson EFT we have 
that while perturbative OPE requires two counterterms 
(Coq and Cot), for non-perturbative OPE this number 
is at least five [sO]. As the shifts in the binding energies 
are only weakly cut-off dependent, we do not require new 
counterterms and we can confidently conclude that OPE 
is perturbative. 

We mention in passing that an alternative possibility 
to check whether OPE is perturbative is the description 
of bound state properties that depend on the existence 
of a _D-wave component of the wave function, which is a 
typical signature of the tensor force. It turns out that D- 
wave probabilities of the molecular states are quite small 
(from 1 — 4%, with a strong cut-off dependence^). Unfor- 
tunately, it looks difficult to find experimental observa- 
tions that could depend on (and consequently constraint) 
the Z?-wave components of the molecular wave- functions. 
Nevertheless, the small size of the I?- wave probabilities 
is consistent with the expectation of them to be a second 
order correction in perturbation theory. 

C. Particle Coupled Channels 

As previously discussed particle coupled channel ef- 
fects are suppressed by two orders in the EFT expan- 



^ One should keep in mind that the D-wave probabiUty is not 
per se an observable quantity, meaning that we should not be 
worried by its moderate cut-off dependence. 
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E-iT/2 (A = 0.5 GeV) 


E-iT/2 (A = 1 GeV) 


Exp f7] 


0++ 


DD,D*D* 


3658 


3669 




1++ 
1+- 


D*D 
D*D,D*D* 


Input 
3730 


Input 
3739 


3872 


0++ 

1+- 

2++ 


DD,D*D* 
0*D* 


3917 - 1 23 
3979 - 1 24 
4012 


3917 - i 50 
3979 - 1 39 
4012 


3917 ± 3 - 1 28tg° 
3942 ± 9 - f 37t 1^ 



TABLE III. Predicted masses and widths (in MeV) of the X(3872) HQSS partners when coupled channels effects are included. 
The contact terms are adjusted to reproduce the X(3872) and X(3915) masses, while OPE effects are neglected. We find 
Coa = -4.16 fm^ and Cot = 2.21 fm^ and Coa = -1.14 fm^ and Cob = 0.35 fm^ for A = 0.5 and 1 GeV, respectively. 





HH 


E - ir/2 (A = 0.5 GeV) 


£;-ir/2 (A = 1 GeV) 


Exp \7] 


0++ 


OD,0*D* 


3690 


3694 




1++ 
i+- 


0*D 
0*D,0*D* 


Input 
3782 


Input 
3782 


3872 


0++ 

1+- 

2++ 


OD,0*D* 
D*D,D*D* 
D*D* 


3939 - 1 12 
3984 - 1 17 
4012 


3937 - 1 31 
3982 - 1 29 
4012 


3917 ± 3 - 1 281:9° 
3942 ±9-5 37^17 



TABLE IV. Predicted masses and widths (in MeV) of the (3872) HQSS partners when coupled channels effects are included. 
The contact terms are fixed to the values given in the caption of Table |I] (i.e., they are adjusted to reproduce the X(3872) and 
X(3915) masses neglecting coupled channel effects). Moreover, OPE interactions are not taken into account either. 

sion. From this we expect coupled channel effects in the 
binding energies of the 0++ and 1^ states to scale as 

\AEB\^\Es\(^]^y , (35) 

where i^s is the binding energy, 75 = ^/—2^Eb the wave 
number of the bound state and Ac the typical momentum 
scale of the coupled channel under consideration, which 
we consider to be a hard scale Ac ~ Ag. The estimation 
above translates into an uncertainty of around 30 MeV 
(40 MeV) for the 0^+ (1^ ) coupled channel, where we 
have employed the wave number of the deepest bound 
state within the coupled channel. It should be noticed 
that in the case of the D*D* molecular states the energy 
shift is complex, as the 0++ (!"' ) state can decay into 
a DD (D*D) meson-antimeson pair. 

In contrast to the OPE corrections, the counterterm 
structure stemming from HQSS is not expected to be 
able to absorb the kind of divergences associated with the 
coupled channel calculations, meaning that the actual er- 
ror in the calculation will probably saturate the previous 
bound. We can check the EFT a priori estimates given 
above by means of a concrete calculation in which the 
particle coupled channel effects are fully taken into ac- 
count. However, as we will see, this task is not trivial, 
specifically in what regards to the choice of the appro- 
priate regularization scheme. To illustrate this point we 
can study the perturbative estimate of the binding en- 
ergy shift induced by coupled channel dynamics. If we 



consider a small change in the potential 
V -^V + 5V , 



(36) 



the perturbative correction to the binding energy is ex- 
pected to be 



5E = {^5V\'^) 



(37) 



where Y^) is the wave function of the bound state. In the 
case of particle coupled channels, the 5V operator reads 



SV^^Vo,fiGQ,p{E)Vp^, 



(38) 



where a represents the channel we are interested in, j3 ^ 
a the other channel and Vap the transition potential from 
channel a to /3, which is proportional to the C06 contact 
operator. We can distinguish two cases, depending on 
whether the unperturbed energy E is (a) above or (b) 
below the j3 channel threshold. The most interesting case 
is (a), corresponding to the modification of the O"*""*" or 
1^ D*D* molecular state energies {a channel) by the 
0++ DD or l+~ DD* states (/3 channel), which lie in the 
continuum for the energies relevant for the a— channel 
states. A direct calculation yields 



6E^ 



Ca/^(A) 

A (27r)3 H 



A (2^ 
2^ 



n 2 



*a(p) 



+ ie 



(39) 
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where ^a{p) is the wave function of the a bound state, 
A:| = —7^ + Aq the momentum of the heavy meson- 
antimeson pair above the threshold, 7q the wave number 
of the a bound state below threshold and Ac the coupled 
channel momentum scale. The integrals are assumed to 
be regularized with a cut-off A and an arbitrary regu- 
lator function that we have not specified yet. The Cap 
transition contact operator can be identified with 2 Cot 
(V^Cob) in the 0++ (1^ ) molecular state. 

As can be seen, if A < fc^, a natural thing to expect 
if the momentum separation of the coupled channels is 

a hard scale, the perturbative correction to the binding 

,2 

/3 



energy is effectively suppressed by a factor of kg ^ Aq 



■/3 



Ca^iA) 



(2^ 



i2n) 



(40) 



that is, by two powers in the counting. However, we can 
also appreciate that the 6Ea correction is strongly scale 
dependent. On the one hand, if we consider that the 
wave function behaves as 



1 



7S 



(41) 



we see that the integral in the first line of the expression 
for the energy shift of Eq. (15^]) diverges as A: 



d^p 
(2^ 



*a(p) oc A. 



(42) 



On the other, the integral related to the decay into the 
continuum state /3 in the second line of Eq. (jJH) di- 
verges as A. However, had we re-expanded the propa- 
gator l/(fc2 _ g-Z) 

in inverse powers of fc^ (as mandated 
by power counting) , the power-law divergence would have 
worsened to A'^, see Eq. PU)) . Putting all the pieces to- 
gether (and ignoring the propagator re-expansion), the 
total divergence in the energy shift is given by C^^ A'^. 

This divergent behaviour tell us that we are required to 
include counterterms to absorb them 0. The problem is 
that the counterterms renormalizing the coupled channel 
dynamics are higher order. In principle the Cob operator 
could do the job, but we need to take into account that 
this contact operator is already determined by the con- 
dition of reproducing the binding energy of a molecular 
state. Thus, we do not expect Cot to balance for the par- 
ticle coupled channel effects. The renormalization group 
behaviour of Co& is approximately given by 



Cofc(A) oc — , 
fiA 



(43) 



It is important to comment at this point that the different diver- 
gences we are discussing correspond to particular choices of how 
to expand in terms of power counting. The full, non-perturbative 
coupled channel calculation is free of divergences. 



for large A, see Eq. (|13l) . This means in turn that Cof, 
can absorb the piece proportional to A^ of the coupled 
channel divergence. 



CapiA) 



d'^p 
(2^ 



(44) 



since Cap is proportional to Cofc. Thus, we are left with 
a residual cut-off dependence of A in the best case. The 
worst case scenario is however when the cut-off and the 
coupled channel scale coincide, A ~ Ac, in which case the 
/3-channel integral peaks. This can be easily appreciated 
if we use a sharp cut-off regulator 



d^q e{A^\q\) 



(2^)3 fc2 _ 



A + i-kpe{A- kp) 



kp 

ylog 



A-kp 



A- 



«/3 



(45) 



where the real part of the integral diverges at A = fc/3 ~ 
Ac, a very puzzling situation (see a related discussion 
in Ref. [5a|)- Of course, the problem can be solved by 
using a smoother regulator, as the gaussian scheme that 
we have been employing along this work. In this case the 
real part integral will show a maximum (but not diverge) 
at A ~ Ac- This signals the transition from a power 
counting in which Ac is a hard scale to a different one in 
which it is a soft scale. 

All this indicates that one should presumably add new 
counterterms at 0{Q) to soften the cutoff dependence 
and make the EFT renormalizable again. At this point 
it is worth mentioning that the EFT treatment of cou- 
pled channel dynamics has been only discussed for the 
case in which Ac is a soft scale [s^, [s^l- However, the 
corresponding analysis for the Ac ~ Ag case has not been 
done yet and will be left for future research [58^ . Indepen- 
dently of the exact form of the power counting for coupled 
channel dynamics, it is clear that higher orders will intro- 
duce new unknown constants that cannot be fixed at the 
moment owing to the scarce experimental data available. 
Nevertheless, here we will present full non-perturbartive 
results including coupled channel effects. Even though 
the energy shifts thus obtained will be cutoff (and regu- 
lator) dependent, they will not vastly deviate from the a 
priori estimates of Eq. psp . reinforcing the (qualitative) 
reliability of the LO predictions. Of course, had we in- 
cluded all the relevant 0{Q) counterterms (and known 
the entire experimental spectrum of D^*-'!)'^*^ states), the 
deviations would have decreased. Furthermore, we no- 
tice that coupled channel effects produce changes in Coa 
and Cob comparable in magnitude to those we should ex- 
pect from violations of HQSS, that is, about 15% for the 
charm quark mass. 

The non-perturbative calculation of the coupled chan- 
nel effects is presented in Table Hill As in previous cases, 
we have adjusted the Coa and Cob counterterms to re- 
produce the A'(3872) and X{3915) masses. We have 
searched for the poles of the scattering amplitude in the 
first and second Riemann sheets (we refer to Ref. [591 
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for further details on this subject). The former are to 
be interpreted as bound states, while the later corre- 
spond to the location of resonant states in the com- 
plex plane, where the real part of the pole position is 
the mass and the imaginary part is half the decay width 
(iJpoio = Af — I r). As can be seen in Table Hill the loca- 
tion of the 1+^ DD* and D*D* states have been shifted 
by about 80 - 85 MeV and 25 MeV respe_ctively. The cor- 
rection to the binding energy of the DD* state is large, 
saturating and even exceeding the EFT expectation. In 
the case of the 0+"'' DD partner of the X(3915), the shift 
in the position of the state is of 40 — 50 MeV, of the 
order of the EFT expectation. Had we use a sharp cut- 
off instead of a gaussian one, the location of the bound 
and resonant states would have drastically changed for 
A = 0.5 GeV. This is not surprising in view of Eq. p5)) 
and the related discussion. However, for larger cut-offs 
such as A = 1 GeV, variations are significantly smaller 
and the results are similar to those obtained in the gaus- 
sian cut-off scheme. 

To further check the uncertainties affecting our results, 
we have also considered the alternative option of using 
the counterterm values of the original uncoupled calcula- 
tion to estimate the coupled channel effects, in which case 
we obtain the results of Table IIVI In this second scheme 
we observe that the change in the position of the states 
agrees much better with the EFT expectations: the ener- 
gies of the 1"' DD* and D*D* states consistently change 
by about 35 MeV, while in the 0++ DD and D*D* states 
we end up with an energy shift of about 20 — 25 MeV. 
The reason for the additional stabilization of the calcu- 
lations may be that we are not forcing the reproduction 
of the Ar(3915) in a cut-off window in which we may 
not expected to obtain this state (owing to the large cou- 
pled channel corrections). In this case the X(3915) state, 
which we do not adjust now, shifts is mass to around 
3940 MeV, with a width of about 15 - 30MeV, depend- 
ing on the value of the cutoff. If we consider that the 
uncertainties coming from the l/mg corrections are of 
the order of 15-30 MeV, the properties of the 0++ D*D* 
state could certainly be accommodated with the existing 
experimental data for this resonance (M — 3917 ± 3 and 
r = 28t^°MeV 0). 



V. HQSS AND DECAY PROPERTIES 

The dynamics of the molecular states studied in this 
work is solely determined, within our approach, by the 
re-interaction of the open charm channels D*^*)!)*^*^. We 
have ignored hidden charm channels like, for example, 
the J/'^uj or ijcUJ. We expect these latter channels to 
have little effect on the inner structure and masses of 
the molecular states, as sugges ted by explicit calcula- 
tions performed in Refs. |60l - l62| . Yet, within the EFT 
approach, it has been customary to ignore hidden charm 
channels in the study of the A (3872) resonance, see e.g. 
Ref. 36] . Analogously, the hidden bottom channels have 



also been ignored in the recent studies [37|, |38j of the 
Z(10610) and Z'(10650) molecular states in the bottom 
sector. Nevertheless, the hidden charm channels can play 
an important role in the decay of some of the states de- 
scribed here, especially if they are placed below the open 
charm D'^*''D(*) thresholds. Moreover, the J/'^ meson 
provides a clear experimental signature and thus its de- 
cay modes are often used in the detection of the XYZ 
states. The detailed study of the hidden charm decays of 
the molecular states described here is beyond the scope 
of this work and we left it for future research. 

However, the generic decay properties of the molecular 
states can be discussed at the qualitative level in the 
basis of HQSS. If we ignore phase space effects, HQSS 
predicts 37, [sl] for the total widths: 



r(i 



= r(2++) 

3 



DD 



(0^ 



rDD-(i" 



2 

rD.D.(i 



■rD.B.(o 



(46) 
(47) 



where we denote each molecular state by its quantum 
number J^^ and additionally its particle content if nec- 
essary. As noticed in Ref. [S^, the relations above can 
also be obtained within the EFT framework we advocate 
by promoting the Cqq and Cot couplings to complex val- 
ues. In this way, one can implicitly take into account the 
multiple decay channels of the molecular states (as with 
an optical potential). In contrast with the bottom sector, 
where the previous relationships were derived, we expect 
however noticeable corrections to Eqs. and (H71) in 
the charm sector. The reason is that both HQSS viola- 
tions and phase space corrections are larger in the charm 
sector than in the bottom one. The relations above in- 
volve total widths and do not necessarily hold for decays 
into open charm channels, where phase space corrections 
are crucial and indeed forbid some decays. For instance, 
if we pay attention to Eq. P5)) . and since rj3]3(0"''+) = 
for an open charm decay into D*D* 0, we will have to 
conclude that r(l++) = r(2++) = rD.s*(0++) = 0. 
However, we find a partial decay width of the order of 
tens of MeV for the A(3915) state into DD. 

Nonetheless, a clear implication of the relationships 
above is that the X(4012) should be a relatively nar- 
row state, just like the A (3872). In addition, if we as- 
sume r(l++) and r(2++) to be much smaller than the 
other decay widths, then we can estimate the total width 
of the A(3710) resonance to be a third of the A(3915) 
width (rD.D.(0++) = 281:1° according to the PDG 0), 
yielding ri3B(0++) - 10 MeV. For the A (3815) and 
A (3955) resonances the situation is similar: HQSS with- 
out phase space considerations predicts them to have the 
same width, but if one takes into account the large D*D 



^ Note that in the infinitely heavy quark limit the D and D* mesons 
are degenerated. Thus, the D*D* decay channel could be open 
depending on the binding energies. 
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contribution to the X(3955), the X{3815) should be nar- 
rower than its partner. 



VI. DISCUSSION AND CONCLUSIONS 

In this work we have argued that the apphcation of 
HQSS to the charmed meson-antimeson system, com- 
bined with the identification of the X(3872) and X(3915) 
resonances as isoscalar DD* and D*D* molecules, im- 
plies the existence of four molecular partners of these two 
states (Table m. This prediction is subjected to a series of 
uncertainties, namely the approximate nature of HQSS 
(especially for the charm sector), the effect of the OPE 
potential and the impact of the particle coupled channel 
dynamics. We have estimated the size of these correc- 
tions within the EFT framework and concluded that the 
HQSS pattern of molecular states is rather stable (Ta- 
bles [TTl [HT] and |IV| . In contrast, the exact location of the 
molecular partners is subjected to moderate uncertain- 
ties of up to 40 — 50 MeV for the most bound cases, in 
agreement with the EFT expectations. 

The determination of the D(*)D(*) family of bound 
states hinges on the assumption that the X(3872) and 
the X(3915) states are molecular. In this regard we find 
it worth commenting that, while the identification of the 
X(3872) as a 1++ loosely bound DD* state is a widely 
accepted hypothesis, the case for the molecular nature 
of the X(3915) is less compelling but nevertheless still 
compatible with the experimental information available 
for this resonance. Thus we expect the conclusions solely 
derived from the X(3872) to be more solid and less spec- 
ulative than those depending on the X(3915). 

In this regard the tentative 2++ D*D* partner of the 
X{3872), which we have called the X(4012) in reference 
to its predicted mass, see Table U is probably the most 
robust and model independent prediction of the present 
work. The X(4012) is not affected by particle coupled 
channel effects and its mass only varies mildly, by about 
2 — 3 MeV, when the OPE potential is included. Perhaps 
in the real world the 1 / mg effects may be larger than we 
have estimated or there may be a further and unexpected 
subleading correction that turns out to be large. In this 
case the X(4012) state might move slightly up above the 
D*D* threshold and become virtual or might descend to 
a lower mass region. Be as it may, we are quite confi- 
dent about the existence of a molecular state with these 
quantum numbers close to the D*D* threshold. 

The prediction of new D(*)d(*) states beyond the 
X(4012) requires the identification of a further XYZ state 
(besides the X(3872)) as a charmed meson-antimeson 
molecule. The X(3915) is a good candidate, which we as- 
sume to be a 0++ D*D* bound state. Of course we notice 
that the molecular interpretation of the X(3915) [T7l - [T9| . 
while plausible, is not so well-established. Consequently 
the three additional 0++ DD and 1+^ DD* and D*D* 
states we obtain from the X(3915), which we call the 
X(3710), X(3820) and X(3955) respectively, see TableU 



should be granted a more conjectural status. Neverthe- 
less, we notice that the only necessary condition for the 
existence of molecular states different than the X(3872) 
and X(4012) is that Cot > 0. 

Other theoretical approaches have also predicted sev- 
eral D*D* molecular-like states, but usually with a mass 
spectrum incompatible with HQSS. In the quark model of 
Ref. (63j there are six hidden charm diquark-antidiquark 
states arranged in a pattern similar to the one we find. 
In particular there is a tetraquark 0++ state at 3723 MeV 
that could be identified with the X{3710) state we ob- 
tain. However, the 2++ state appears at 3952 MeV and 
is identified with the X(3940) resonance [13, HH. Un- 
less a considerable violation of HQSS is taking place, 
this tensor state is too tightly bound to be considered 
the HQSS partner of the X(3872). Curiously, the two 
1^ hidden charm diquark-antidiquark states of Ref. [gl] 
are located at a similar depth below the D*D* and DD* 
thresholds respectively and therefore respect the HQSS 
expectations. 

Another interesting theoretical approach for the study 
of hidden charm resonances is the hidden gauge formal- 
ism, using an extension of the SU (3) chiral lagrangians 
to SU (4) that implements a particular pattern of SU (4) 
flavor symmetry breaking. Within this framework. Gam- 
merman et al. [lO] have obtained a 0++ DD molecular 
state in the vicinity of 3700 MeV, that is to be identified 
with the X{3710) DD molecular state we predict. The 
extension of the hidden gauge to axial states [6l| predicts 
(among others) a negative C-parity state at 3840 MeV, 
not far way from the 3815 — 3820 MeV mass range we 
obtain for the 1^ DD* state. Finally, the related ex- 
ploration of resonances generated by the interaction of 
two vector mesons in Ref. [g^l predicts a series of 0++, 
1+" and 2++ D*D* states. The 0++ D*D+ resonance 
is found in the region around 3940 MeV. Though not 
identical, this figure does not differ much from the mass 
of the X(3915) resonance that we employ as input. The 
1+^ D*D* state of Ref. [12] matches rather well with 
the mass of the X(3955) state we obtain. However, 
the 2++ D*D* isoscalar resonance is considerably dif- 
ferent from the X(4012) state: its mass and width are 
M = 3929 ± 3 MeV and F = 29 ± 10 MeV respectively, 
where the dominant decay channel is DD. This mass, 
which is clearly incompatible with the HQSS pattern, is 
the result of the remarkably strong vector- vector interac- 
tion generated by the hidden gauge model. Curiously, the 
properties of this tensor D*D* resonance are strikingly 
similar to those of the Xc2 {'^P) charmonium state l65| 
(sometimes referred to as the Z(3940)): M = 3927.2±2.6 
and F = 24 ± 6 MeV, decaying mostly to DD 0. 

The comparison of the HQSS spectrum with experi- 
mentally known states is however incomplete. In prin- 
ciple there is so far no experimental evidence in fa- 
vor (or against) of the positive C-parity X(3710) and 
X(4012) states we predict. Interestingly, the properties 
of the predicted X(3955) molecular state are not very 
different from what is experimentally known about the 
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X(3940) resonance [13, [U, which have been observed 
to decay into DD* (just as would have been expected 
for a 1^ D*D* state). There is a problem though in 
this identification: the X(3940) is suspected to be a pos- 
itive C-parity state, while the X{3955) has negative C- 
parity. The reason is that the usual production mecha- 
nism e+e~ — )■ 7* — >■ J/'i/X favors the generation of pos- 
itive C-parity XYZ states, owing to the quantum num- 
bers of the intermediate virtual photon and the final J/4'. 
This mechanism also implies that any prediction about 
negative C-parity states will be more difficult to con- 
firm or discard experimentally. However, even though 
not so probable, the production of the final J/'^X state 
may happen via two virtual photons, in which case the 
XYZ resonance may have negative C-parity. This al- 
ternative production mechanism is not always as sup- 
pressed as expected, as demonstrated in Refs. |52l.[53j for 
e^e~ — ^ 7*7* — > J J/4'. In principle, a similar mech- 
anism could take place in the X(3940) state, in which 
case the identification with a 1^ DD* molecule would 
be very appealing, but it may also be possible that the 
X(3955) molecular state has simply not been observed 
yet. 

We have also examined the role played by the OPE 
potential in the D^*)©^*) system (Table HI)) . In agree- 
ment with the conclusions of Ref. [30(], we have verified 
that pion exchanges can be treated perturbatively in the 
case of isoscalar charm meson- ant imeson molecules. Cu- 
riously, the suppression of the OPE effects is larger than 
naively expected in terms of the power counting. The 
remarkable simplification is that the EFT framework we 
are proposing consists entirely on contact range inter- 
actions at lowest order. According to Ref. |3(| this in- 
teresting simplification will also apply in the isovector 
charm sector, as well as the isovector bottom one [38l |. 
In this regard non-perturbative OPE seems to be only 
required in the case of isoscalar bottom meson- antimeson 
molecules. This possibility, and the corresponding EFT, 
was partially explored in Ref. [11] for the BB*/ B*B case. 
Lastly, particle coupled channel dynamics are suppressed 
by two orders in the counting, as expected, but a more 
complete analysis would be welcomed, specially in what 
regards to regulator dependence. 

Even though all the previous states and their possible 
identification with theoretically predicted / experimen- 
tally known resonances are contingent on the validity of 
the molecular hypothesis for the X(3915), the bottom- 
line of the approach we advocate is that, provided we 
identify at least two molecular states representing two 
different combinations of the Coa and Cob counterterms, 
we will be able to predict the full molecular spectrum of 
theD(*)D(*) system. IftheX(3915) identification proves 
erroneous in the future, the finding of a different molec- 
ular state candidate could be used to obtain the remain- 
ing states. Owing to the contact range character of the 
present EFT framework at lowest order, the calculational 
effort involved in this task will be minimal. 



ACKNOWLEDGMENTS 

We thank E. Ruiz Arriola for discussions. This re- 
search was supported by DGI and FEDER funds, un- 
der contract FIS2011-28853-C02-02, and the Spanish 
Consohder-Ingenio 2010 Programme CPAN (CSD2007- 
00042), by Generalitat Valenciana under contract 
PROMETEO/20090090 and by the EU HadronPhysics2 
project, grant agreement no. 227431. 

Appendix A: The Effective Lagrangian at Lowest 
Order 

In this appendix we write the EFT Lagrangian that 
describes the strong interactions of heavy mesons and 
antimesons containing a heavy quark Q or antiquark Q 
respectively. We use the matrix field iJ^'^) (iJ^'S)) to 
denote the following combination of the pseudoscalar and 
vector heavy-meson (antimeson) fields 

HP - ^ {p:fh' - pPi.) , (Al) 

HiQ)- ^ (p;(0)a^M _ p(Q)-^5) i_i , (A2) 

where the pseudoscalar meson (antimeson) fields are rep- 
resented by pP {PP), while P*^^^ {Pa^^^) is employed 
for their vector HQSS partners (see, for example, Ref. [6g| 
for further details). Finally, v is the velocity parameter. 
In principle there should be a w subscript to indicate that 
we are defining the fields for a specific value of u, but we 
have omitted it to avoid complicating the notation. The 
fields are isospin doublets (hence the index a), where for 
the pseudoscalar meson and antimesons we have 

P^Q/Q) ^ (p0^p+)^ (A3) 
pW/Q) ^ (po^p-)^ (A4) 

plus the analogous expressions for the vector case. The 
heavy quark/antiquark superindex changes depending on 
whether we are considering charm or bottom meson fields 
(in the charm case, D'^ and contain the quark field, 
while in the bottom case B*^"^ and i?+ contain the anti- 
quark field) . The heavy vector meson and antimeson are 
subjected to the additional condition 

V ■ P<^^ = , (A5) 
V ■ P*^^^" = , (A6) 

which in turn defines the three different polarizations of 
the heavy vector mesons. 

The fields 7?!'^^ and iJ^'^)" respectively transform as 
a (2, 2) and (2, 2) representation under the heavy quark 
spin (g) SU(2)y isospin symmetry [gT], that is 

iJ^^) S'(ij('5)j7l') , (A7) 

H^Qy ^ {uh'^^^Y S\ (AS) 
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where S is the heavy quark spin transformation and U 
the isospin one. The hermitian conjugate fields are 



and transform as [63| 



(A9) 
(AlO) 

(All) 

(A12) 
I 



Of course, the Lagrangian should be invariant under the 
previous symmetry transformations. 

At leading order in the EFT expansion the Lagrangian 
can be written as the sum of two contributions 



(0) 
4H 



C. 



(0) 

ttHH 



(A13) 



where the first one contains a 4-meson interaction vertex 
and the second the meson-pion vertex. The 4-meson con- 
tact range Lagrangian consistent with HQSS and chiral 
symmetry 36] reads: 



^fn ^ -^Oa Tr 
+ Sn„ Tr 



Tr 



Tr 



i^OfcTr 



Tr 



7m 75 



Tr 



Tl5 



(A14) 



where Tab are the Pauh matrices, and a,6,r and s are 
isospin indices. We notice that for each isospin channel 
(/ = 0, 1) we have only two independent constants. 

On the other hand, at leading order in the chiral ex- 
pansion the HHtt and HHtt couplings are determined by 
the Lagrangian [g^l 



(0) 

ttHH 



V2U 
+Tr 



Tr 



^(Q)b^(Q) 



Tl5 



-0(^2) 



(A15) 



where tt is the relativistic field that describes the pion, 
g is the PP*tt coupling and ~ 132 MeV the pion de- 
cay constant. In the charm sector, g has been deter- 
mined from the D* meson decays in Refs. [i^ |4l| yield- 
ing g — 0.59 ± 0.01 ± 0.07, which we approximate by 
(7 ~ 0.6. In the strict heavy quark limit, the latest lattice 
QCD resuhs suggest the value g = 0.449 ±0.047 ±0.019, 
see Refs. [1^ [69| . In the normalization above the pion 
field has dimensions of [energy], while the heavy meson 
or antimeson fields ff^'^) or H'^'^^ have dimensions of 
[energy]^^^: as usual in heavy quark physics, we employ a 
non-relativistic normalization for the heavy mesons that 
differs from the usual relativistic convention by a factor 
of y/Mn (see for instance Ref.fs^). 



Appendix B: Projecting the Potential into the 
Partial Wave Basis 



In this appendix we delineate how to project the heavy 
meson-antimeson potential into the partial wave basis. In 
first place we define the non-relativistic potential for the 
transition (not necessarily elastic) 

H(1)H(2) ^H(1')II(2'), (Bl) 



in terms of the tree level scattering amplitude 

7;roo = -*V(l + 2^l' + 2'), (B2) 

where 1,2 and l',2' schematically represent the initial and 
final state of each of the particles. For heavy meson- 
antimeson scattering the initial (final) state is completely 
represented by the momentum exchanged between the 
particles p = Pi — P2 (p ' = Pi' — P2') and by the total 
and third component of the spin of each of the parti- 
cles, which we can collectively call a = {{Simi){S2m2)} 
{a' = {{S[m[){S^m'2)}). If we compute V(l + 2 ^ l'+2') 
in terms of the usual Feynman rules (each vertex con- 
tributes with iC, additional i factors for each pion prop- 
agator in the case of OPE contributions, etc.), the re- 
lationship between the non-relativistic potential and the 
invariant scattering amplitude is 

{p';a'\V\p;a) = ^V{p,a ^ p' ,a') . (B3) 

Notice that the factor dividing the invariant scattering 
amplitude is 4, instead of the usual 4^/M^JhM{Ml, 
owing to the \/Mh normalization factor included in the 
heavy meson/antimeson fields. 

Now we specify the procedure for the partial wave pro- 
jection of the potential. To denote the different par- 
tial waves we employ the spectroscopic notation 
where 5, L and J are the intrinsic, orbital and total an- 
gular momentum. With this, we define the states with 
good angular momentum as follows 



1 



J2 {LSJ\MlMsM) 



Ml Ms 



\p;JMLS) 



where p is the modulus of the center-of-mass (cm.) mo- 
mentum p (= p X p) of the HH pair and {LSJ\MlMsM) 



X / dn{p)YL,MApMSMs), (B4) 
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is a Clebsch-Gordan coefficient. Tfie normalization of tlie 
states above can be determined from the normalization 
of the plane wave basis, that is 

(p"; S'M's \p; SMs) = {2tt)' S'{p- p') Ss,s'SmsM's > 

(B5) 

yielding 

{p'-J'M'L'S' \p;JMLS) = 

(2^)3 y Sjj,Smm'Sll'Sss' ■ (B6) 
I 



In this basis the partial wave projection of the potential 
reads 



vfAip'.P) = {P; JML'S' \V\p ; JMLS) 

= 1. j dn{p) j dn{p') E {LSJ\MaW{L'S'J\M'^M'sM)Yl. j,,,^{p')YLMM 

MlMsM'^M'^ 

X J2 (^i^2^|mim2Ms)(5i5^5'|mW2A^5)(p';(^W'i)(^>2)l"t^|p;(^i^'^i)(^2m2)) , (B7) 

mim2m'j^m2 



where thanks to rotational invariance the above matrix 
element is independent of the third component of the 
total angular momentum M. 

Appendix C: The Lowest Order Heavy 
Meson- Ant imeson Potential and Its Partial Wave 
Projection 

The lowest order HH potential contains a contact and 
a finite range (pion exchange) piece. The EFT poten- 
tial can be derived from the tree level scattering ampli- 
tudes resulting from the C^^^ and C^^j^j^ Lagrangians of 
Eqs. (jA14p and (jAlsp . Even though the partial wave 
projection of the contact piece is trivial, we will start 
with the OPE potential in order to fix the notation. The 
tree level invariant amplitude can be obtained from the 
pion-meson dfjju Lagrangian of Eq. (jA15[) . taking the 
schematic form in the strict heavy quark limit, 

VoPEb-p)oc ^l-ff -p , (CI) 

where q = p — p' (that is, the potential is local) and a, 
h is the corresponding polarization operator in each of 
the HHtt vertices. The proportionality factor is /^f"^ 
times a sign that depends the pseudoscalar or vector 
nature of each of the particles in the initial and final 
states. The modifications to take into account in the 
above equation the mass difference between the pseu- 
doscalar and vector heavy mesons masses are also dis- 
cussed in Refs. [13, H^l- Since we will not be considering 
particle coupled channels with the OPE interaction, this 
becomes an issue in this work only for the DD* — )■ D*D 
channel. In that case in Eq. (jCll) . should be sub- 



stituted by /x^ = — (m£). — mo)^. Indeed, fi^ < 0, 
since there is a very small absorptive contribution from 
the DDtt channel. We will neglect it, as in Ref. [2^ . 
and in that case we will consider only the real part of 
the potential. We now continue by Fourier transforming 
the amplitude above into coordinate space. This step, 
though counter-intuitive at first sight, will enormously 
facilitate the calculation of the partial wave projection of 
the potential in momentum space. We remind that the 
coordinate and momentum space potentials are related 
by 

{p';a'\V\p;a) = J d^r e^^P-P'^-Ua'\V{r)\<j) , (C2) 

and we use the symbols a and a' to encode all the spin 
indices (see the previous appendix). Then we make use 
of a well-known relationship 




from which we obtain the form of the OPE potential 

VopE{r) (X (5^(f) - (vc{r) d-b + vrir) 5'i2(a, b)) . 

(C4) 

In this equation Si2{d,b) is the tensor operator, which 
we define as 

J. 3(a ■r)(b-r) ^ - 
bi2[a,b) = ^ a • 0. (05) 
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In turn, the central and tensor pieces of the potential, Vc 
and vt, are given by 



vc{r) = — 
LZtt 



vrir) = vcir) 1 + 



(C6) 
(C7) 



Of course, we are interested in the partial wave projection 
of the potential above, Vjj^fi^{r), since its partial wave 
Fourier transform provides the multipole expansion of 
the potential in momentum space (Eq. (jB7|) ) 



L-L' 



+00 



dr r^jL [p r) jl' {p'r) Vfi?L [r) 



(C8) 

where jhix) represent the spherical Bessel function of 
order L. The advantage of this expression is that it can 
be analytically evaluated with relative ease for the OPE 
potential. After a bit of Racah algebra and taking into 
account all the signs and factors we have obviated so far 
(again, details can be consulted in Ref. [lOl), we arrive at 
the final expression for the OPE potential in the partial 
wave basis 



{VoPEfjldT) = 



Cr 



9^ S{r) 



(C9) 



where all the calculational complications are conveniently 
hidden in C12 and S'12, the partial wave projections of the 
a • h and 5*12(0,6) operators: C12 and 5*12 depend on J, 
i, L' , S and thus encode all the information required to 
determine the coordinate space potential in a particular 
partial wave. 

The most compact way to write the C12 and 5*12 fac- 
tors is in matrix form, where the matrix is defined in a 
basis formed by the set of partial waves with well-defined 
quantum numbers J^*^. We can illustrate this by con- 
sidering all the J^'-' combinations that contain S-waves, 
that is, the J^'~^ values we have studied in this work. 
First we consider the set of 0+"'" partial waves, defined as 



6(0++) = [DDCSo),D*D*CSa),D*D*{'>Da)^ , 



from which the C12 and S12 matrices are 



Ci2(0 



Sl2(0^ 




(CIO) 



(Cll) 



(C12) 



where we have only shown the lower and diagonal com- 
ponents as the matrices are symmetric. Next we move to 
the 1++ case, for which we have the particle states DD* 
and D*D that we need to arrange in states with good 
C-parity. We thus define 



[DD*{Tf)] = -= [DD* -riD*D] . (C13) 
v2 



In this convention the intrinsic C— parity of these states 
is independent of the isospin and equal to 77. The 1++ 
basis thus reads 



S(l++) = {[DD*{+)]CS,), [DD*i+)]CD,), 

D*D*{^Di)}, (CM) 

for which we obtain the matrices 



(C15) 



(C16) 




Sl2(l++) 



The next case is 1+ , in which we employ the basis 

B{1+-) = {[DD*{^)]eS,), [DD*{-)]eD,), 

D*D*{^Si),D*D*{^Di)Y (C17) 



and get the matrices 



Cl2(l + -) 



Sl2(l 



1 
-2 1 
V -2 1/ 



/ \ 

-V2 1 

-\/2 
V-V2 1 -V2 ij 



(C18) 



(C19) 
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The most complex case is 2++, in which OPE mixes ah the possible particle channels. We thus work in the basis 
B{2++) = i^DDCD2), [DD*{+)]{^D2), D* D* D2) , D* D* S2) , D* D* {'^02), D* D*{^G2)} , (C20) 
and obtain the matrices 

\ 



Cl2(2^ 



Si2(2^ 



/ 

-1 

-V^ 2 

-1 

0-1 

V 

/ 



1 









3_ 

35 



3 

12 
V35 



2_ 

^/7 
6 
V35 



-1 / 



(C21) 



(C22) 



12 10 , 

7V5 7 / 



Now we consider the contact range piece of the poten- 
tial, which can be easily derived from the 4-meson C^^^ 
Lagrangian of Eq. (jA14l) . For that, we expand the C^^^ 
Lagrangian into its explicit representation in terms of 
pseudoscalar and vector heavy meson fields (detailed ex- 
pressions can be consulted in Ref. [13] )■ The subsequent 
partial wave projection is straightforward (owing to the 
simplification of working with S- waves only): 

s' s ^(^) 

(^c)jL'Li^) = i^Oa - Eoa fl ■ f2) Sjs Sj'S' Slo Sl'O 

+ ^(Dob-EobTi-T2)Ci2, (C23) 

from which the potentials of Eqs. P^ - (PT|) are derived. 
Notice however that in Eqs. (|T8l) ~ (|2T|) we have speci- 
fied the potential according to the J^^ quantum num- 
ber of the heavy meson-antimeson system. The relation 
between the couplings of Eqs. P^ - (PT|) and the corre- 
sponding ones in the cf^ Lagrangian is provided by the 



expressions 

Coa — Doa + 3 Eoa , 

Cob = Dob + iEob (C24) 

where we have isolated the isoscalar contribution 0. If 
we are considering the full potential resulting from the 
sum of the contact range and OPE contribution 

= iVcfj'Lir) + {VoPEfjldr) , (C25) 
we see that the 5{r) contribution within the OPE poten- 
tial can be absorbed within the contact range piece by 
means of the replacement 

Cob — > Cob — 7775- , (C26) 

which in fact is able to account for the bulk of the change 
of this operator when we add the OPE potential, see the 
related discussion in Sect. IIVI 
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